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Abstract 

The dynamic melting of vortex lattices in type II superconductors is consid- 
ered. A field-theoretic formulation of the pinning problem allows the average 
over the quenched disorder to be performed exactly. A self-consistent theory 
is constructed using a functional method for the effective action, allowing a 
determination of the pinning force and the vortex fluctuations. The phase 
diagram for the dynamic melting transition is determined numerically. In 
contrast to perturbation theory, the self-consistent theory is in quantitative 
agreement with the prediction of a recent phenomenological theory and sim- 
ulations and experimental data. 
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The advent of high temperature superconductors has stimulated a renewed interest in 
the dynamics of vortices in type II superconductors. In this letter we consider the influence 
of quenched disorder on the dynamic melting of a vortex lattice. This non-equilibrium 
phase transition has been studied experimentally as well as through numerical simulation 
and a phenomenological theory and perturbation theory - [^]. The notion of dynamic 
melting refers to the melting of a moving vortex lattice where in addition to the thermal 
fluctuations, fluctuations in vortex positions are induced by the disorder. A temperature- 
dependent critical velocity distinguishes a transition between a phase where the vortices 
form a moving lattice, the solid phase, and a vortex liquid phase. 

We consider a two-dimensional system (normal to n) since we have a thin superconduct- 
ing film in mind, or a 3D layered superconductor with uncorrelated disorder between the 
layers. The description of the vortex dynamics will be based on the Langevin equation 0, 

i 

r/URi + J2 = F - VV(R + unt) + Cm, (1) 

R' 

where is the displacement at time t of the vortex which initially has equilibrium position 
R, and r] is the friction coefficient (per unit length of the vortex). The dynamic matrix, 
$RFt/, of the triangular Abrikosov vortex lattice describes the harmonic interaction between 
the vortices, and is specified within the continuum theory of elastic media by the compression 
modulus, cii, and the shear modulus, cge, 

^ ^ 00 / cnql + CQ^ql (cn - Cee)?^?^ \ 
B \ (cii - CGe)qxqy c^&ql + Cng^ / ' 

where (^q = h/2e is the fiux quantum, and B the magnitude of the external magnetic 
field, B = Bh, and (po/B is therefore equal to the of the unit cell of the vortex 

lattice. The force (per unit length) on the right hand side of eq. (|l]) consists of the Lorentz 
force, F = 0oj X n, due to the (assumed constant) transport current density j, and the 
thermal white noise stochastic force, ^j^^, is specified according to the fiuctuation-dissipation 
theorem (^m^mt') ~ ^rjTSit — t')^oa'^RR'> V is the pinning potential due to quenched 
disorder. The pinning is described by a Gaussian distributed stochastic potential with 
zero mean, and thus characterized by its correlation function {V{x)V{x')) = i/(x — x') = 
z/q/ (27rrp) exp{ — |x — x'p/ (2rp)}, taken to be a Gaussian function with range and strength 
uq in our numerical calculations. 

The average vortex motion is conveniently described by reformulating the stochastic 
problem in terms of the dynamic field theory |^. The probability functional for a realization 
{urjIr of the motion of the vortex lattice is expressed, using the equation of motion, as a 
functional integral over a set of auxiliary variables {uRt}R, and we are led to consider the 
generating functional 

Z[F,3] = fUVnnt [HVhn't' e^'^l"'"], (3) 

R R' 

where in the action, S[u, u] = u{D~^^u + F — W + C) + Ju, matrix notation is used in order 
to suppress the integrations over time and summations over vortex positions and Cartesian 
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indices. The retarded Green's operator is given by — D^^u = r/URt + X]r' ^RR'UR't, and its 
Fourier transform is the matrix in Cartesian space D]^\q. u) = irjujl — $q. The average with 
respect to thermal noise and quenched disorder is immediately performed and we obtain the 
averaged generating functional 



I 



V(f)e 



,iS[<l>]+if<l>+^<t>K(l, 



(4) 



where in order to obtain a self-consistent equation for the two-point Green's function, we 
have added a two-particle source term K. We have introduced the notation = (u, u) and 
/ = (F, J), and the source term, J, coupling to the vortex positions, u, allow us to generate 
the vortex correlation functions, say. 



(URtUR/J/)) = - 



6^Z 



(5) 



K=0,3=0 



The action upon averaging, S = Sq + Sy, consists of a quadratic term, 5*0 [0] = (f)D~^(f)/2, 
specified in terms of the free inverse symmetric matrix Green's function (where the matrix 
in the dynamical, or Keldysh, indices in addition is a matrix in Cartesian indices, and 
time and vortex positions) 



D 



_i ^ ^ 2t7jT6a(s6nR'S{t - t') 



and a term originating from the disorder 




(6) 



(7) 



RR' 



In order to obtain an equation for the pinning force, we consider the effective action 
r[0,G] = W[f,K] - f'4>- tK4>/2 - iTr{GK)/2, the generator of two-particle irreducible 
Green's functions, i.e., the Legendre transform of the generator of connected Green's func- 
tions, W[f,K] — —ilnZ[f,K], which satisfies the equations 



5(f) SG 2 ' 



(8) 



where (p-^^ is the average field, and G is the full connected two-point matrix Green's function 
of the theory, and Tr denotes the trace over all variables. 

In the physical problem of interest the sources K and J vanish, and the full matrix Green's 
function has, due to the normalization of the generating functional, Z[F, J = 0, ^=0] = 1, 
the structure in Keldysh space 



G 



G^ 

QR QK 



((5m" 5uf^)) \ 
{{Su" 5u^)) {{Su'' 5u^)) ) ' 



(9) 



where 5n^t — Ur* — ((urj)) and ^Urj = Urj — ((ur*))- The retarded Green's function G^^ 
gives the linear response to the force F^, and is the correlation function (both matrices 
in Cartesian indices as indicated). 
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As shown in ||^, the effective action can be written on the form r[0, G] = S[(j)] + 
iTiiDs^G - \n{D-^G) - l)/2 - z ln(e^^'"*[*''^l)G ^ where D^^ = 5^S[^]/5^S^, and ^int[0,^] 
is the part of S[(j) + ip] which is higher than second order in ip in the expansion around 
the average field 0. The superscript "2PI" on the last term indicates that only two-particle 
irreducible vacuum diagrams should be kept in the interaction part of the effective ac- 
tion, and the subscript that propagator lines represent G, i.e., the brackets with subscript 
G denote the average (e'^'-^^^'^P^) g = (det iG')-^/^ /D^ giVG-V/2gi5i„t[0,^]_ or^er to get 
a closed expression for the self-energy in terms of the two-point Green's function an ap- 
proximation is called for, and we expand the exponential and keep only the first term, 
— z ln(e*'^'"*f'^''^l)^^^ ^ {Sv[(j) + i-e., the Hartree approximation where diagrams with 

propagators connecting different impurity correlators are neglected (which can also be ex- 
pressed as a Gaussian fluctuation corrected saddle-point approximation 0). 

For the problem of interest, the two-particle source, K, vanishes, and the second equation 
in (H) yields the Dyson equation, G^^ = — S[0, G], with the Keldysh matrix self-energy 
given by 



~ I 5G 



K = ■ (10) 

((u» = 



The Dyson equation and eq. ([T0|) constitute a set of self-consistent equations for the Green's 
functions and the self-energies. Since the expectation value of the auxiliary field vanishes due 



to the normalization of the generating functional, the average field entering eq. (10) is = 
(((uRf)), ((uRt))) = (0, vt), where v is the average lattice velocity. The matrix self-energy has 
two independent components, say and T,^, since S^^(Rt,R't') = Sf^(R't', Rt). From 
eq. (0) we obtain for the case of N vortices Ef^(Rt,Rr) = -i/{Na'^) J2^^u{k)k^kf3e'^^, 
and S-^(q, u;) = (j^(q, a;) — cT^(q = 0, ci; = 0), where the Fourier transform has the Cartesian 
components af^(Rt,Rr) = l/(A^a2) i/(k)A;„A;^(kG^(Rt, Rr)k)e*^''. The influence of 
thermal and disorder induced fluctuations is described by the complex phase = ikMk + 
k ■ (R - R' + v(t - t')), specified by the Cartesian matrix M„^(Rt, R't') = i{G^piRt, Rt) - 
G^p(R.t, R't')). Using the Langevin equation and the first equation in (^ we obtain for the 
pinning force, Fp = —{{VV)), 

oo 

Fp = tJ2 [dt[j^kiy{k){kG''{Rt,R't')k)e'^\ (11) 



R' 



Let us recall the argument for determining the phase diagram for dynamic melting of a 
vortex lattice presented in ^j. There the disorder induced fluctuations were estimated con- 
sidering the correlation function, ^/^^(x, t) = {{fp{yi,t)fa{0,0))), of the pinning force density 
of the vortices, f(x, t) = — ^^(5(x — R^{t))'VV{'x — vt). Neglecting the interdependence 
of the fluctuations of the vortex positions and the fluctuations in the disorder potential, 
the pinning force correlation function factorizes, and since in the fluidlike phase the motion 
of different vortices are "incoherent", one gets /€/3q(x, t) ^ — n^5(x)V/3VQ,i^(vt), being 
the density of vortices. In analogy with the noise correlator the effect of disorder induced 
fluctuations is represented by a "shaking temperature" 



^sh = -; y^mxmt i^aai^, t) = — ^= , (12) 
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where in the last equahty it is assumed that the pinning force is small compared to the 
friction force, i.e., riv ~ F. An effective temperature is then obtained by adding the "shaking 
temperature" due to disorder to the temperature, and according to eq. (|T2|) the effective 
temperature decreases with increasing external force (i.e., with increasing average velocity 
of the vortices). As the external force is increased the fluid thus freezes into a lattice. 
The value of the external force for which the moving lattice melts, the transition force 
Ft, is in 10] taken to be the value for which the effective temperature equals the melting 
temperature in the absence of disorder, Tes{F=Ft) = T^, and has, according to eq. (f2[), 
the temperature dependence Ft = uq/ {4:^/27rrp{Tm — T)), as the temperature approaches the 
melting temperature of the ideal lattice from below (we note that the transition force for 
strong enough disorder exceeds the critical force for which the lattice is pinned Ft > Fc = 

We now describe the calculation of the physical quantities of interest using the self- 
consistent theory. The conventional way of determining the melting transition is to use the 
Lindemann criterion, which states that a lattice melts when the displacement fluctuations 
reach a critical value ((u^)) = c|^a^, where cl is the Lindemann parameter which is typi- 
cally ranging in the interval from 0.1 to 0.2. In two dimensions the position fluctuations 
of a vortex diverge even for a clean system, and the Lindemann criterion implies that a 
two-dimensional vortex lattice is always unstable against thermal fluctuations. However a 
quasi long-range translational order persists up to a certain melting temperature 0. As a 
criterion for the loss of long-range translational order a modified Lindemann criterion in- 
volving the relative vortex fluctuations, (((u(R+ ao,t) — u(R, t))^)) = 2c|a^, where ao is a 
primitive lattice vector, has successfully been employed P], and its validity verified within a 



variational treatment [0. The relative displacement fluctuations are specified in terms of 
the correlation function, (((u(R + ao, t) - u(R, t))^)) = 2itr(G^(0,0) - G'^(ao,0)), where tr 
denotes the trace with respect to the Cartesian indices. The correlation function is according 
to the Dyson equation = G'q^(S^ — 2iriT)G-^^, where the influence of the quenched 
disorder appears explicitly through S^' and implicitly through and S"^ in the advanced 
and retarded response functions. Furthermore, the self-energies depend self-consistently on 
the response and correlation functions. We have calculated numerically the Green's func- 
tions and self-energies and thereby the vortex fluctuations for a vortex lattice of size 8x8, 
and evaluated the pinning force from eq. (pH]). 

We determine the phase diagram for dynamic melting of the vortex lattice by calculating 
the relative displacement fluctuations for a set of velocities, and interpolate to find the 
transition velocity, Vt, i.e., the value of the velocity at which the fluctuations fulfill the 
modified Lindemann criterion. The transition force is given by the averaged equation of 
motion. Ft = rjvt + Fp{vt), and is then determined by using the numerically calculated 
pinning force. Repeating this procedure for various temperatures determines the melting 
curve, i.e., the temperature dependence of the transition force, Ft{T), separating two phases 
in the FT-plane: a high velocity phase where the vortices form a moving solid when the 
external force exceeds the transition force F > Ft{T), and a liquid phase for forces below 
the transition force. 

In order to be able to compare the results of the self-consistent theory to the simulation 
results in , we choose the same parameters. There the melting temperature in the absence 
of disorder is taken to be = 0.007 (2(0o/47rA)^ is chosen as the unit of energy per unit 
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length, where A is the London penetration depth), the value obtained by simulations of 
clean systems, and assumed equal to the Kosterlitz-Thouless temperature, Tkt = Cgeo /47r, 
determining thereby the shear modulus to have the value Cge = 0.088 (a is chosen as the 
unit of length) |Tl|. The London penetration depth is according to @] taken to equal the 
range of the vortex interaction potential used in [Q], which is approximately equal to the 
lattice spacing, ao, giving for the compression modulus cn = IGTrA^cge/cto — SOcee — 4.4 
1^. The range and strength of the disorder correlator in the simulations are Vp = 0.2 and 
uq = 1.42 ■ 10~^ in the chosen units. 

Once the Lindemann parameter is determined, our numerical results for the relative 
displacement fluctuations can be used to obtain the dynamic phase diagram. In order to do 
so we calculate "melting curves," as described above, using the self-consistent theory for a 
set of different values of the Lindemann parameter. We find that these curves have the same 
shape, close to the melting temperature, as the melting curve obtained from the shaking 
theory, T = Ci — C2/Ft. The curve which intersects closest to the melting temperature 
Tm = 0.007, the one depicted in the inset in the figure, then determines the Lindemann 
parameter to be given by the value c\ = 0.0153. The corresponding phase diagram obtained 
from the self-consistent theory is shown in the figure, as well as the simulation data in 0, 
and the melting curve obtained from the shaking theory. The melting of the vortex lattice 
was in the simulations indicated by an abrupt increase in the structural disorder, a different 
melting criterion, and the agreement of the self-consistent theory with the simulation data 
are therefore independently validating the use of the modified Lindemann criterion. The 
shaking theory is seen to be in remarkably good agreement with the self-consistent theory 
even at lower velocities where the shaking argument is not a priori valid, a feature which, 
however, is less pronounced for stronger disorder. It is also of interest to recall, that while 
the melting curve obtained from the shaking theory was based on arguments only valid in 
the liquid phase, i.e., freezing of the vortex liquid was considered, the melting curve we 
obtained from the self-consistent theory is calculated assuming a lattice, i.e., we consider 
melting of the moving lattice. As apparent from the inset in the figure, the critical exponent 
obtained from the self-consistent theory, 1.0, is in excellent agreement with the prediction of 
the shaking theory, where the critical exponent equals one. Furthermore, we find that the 
self-consistent theory yields the value 1.65 ■ 10~^ for the constant C2, a value fairly close to 
the one predicted by the shaking theory, z/o/(4v^27rrp) = 1.77 ■ 10~^. 

The melting curve predicted by lowest order perturbation theory is also shown in the 
figure. The correlation function is in this case given by GqS^^ = -Dq^[Sqi^'' — 2iriT — 

2z?7T(Eqi,^''Dq^+Dq^Sqi,^'')]Dq^, whcrc the self-energies are calculated to first order in z/q. As 
to be expected, the perturbation theory results are at high velocities in good agreement with 
the self-consistent theory. However, we observe that the melting curve obtained from lowest 
order perturbation theory deviates markedly at intermediate velocities from the prediction 
of the self-consistent theory, and thereby the shaking theory, which is known to account well 
for the measured melting curve, Hellerqvist et. al. [IH. 

In conclusion, we have developed a self-consistent theory of the dynamic melting tran- 
sition of a vortex lattice, enabling us to determine numerically the melting curve directly 
from the dynamics of the system. The self-consistent theory corroborates the phase diagram 
obtained by the phenomenological shaking theory, whereas lowest order perturbation theory 
does not. The melting curve obtained from the self-consistent theory is found to be in good 
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quantitative agreement with the phenomenological theory as well as with simulations and 
experimental data. 
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FIG. 1. The dynamic melting phase diagram. The points on the melting curve separates the 
two phases - for forces stronger than the transition force the moving vortices form a solid, and for 
weaker forces a liquid. The dots in the boxes represent points on the melting curve obtained from 
the self-consistent theory for an 8 x 8 lattice, while the three stars represent the simulation results 
in [7]. The crosses represent the results of lowest order perturbation theory. The dashed line is 
the curve Ft{T) = 1.77 • 10-''/(0.007 - T), the melting curve obtained from the shaking theory. 
Inset: Temperature as a function of the inverse transition force obtained from the self-consistent 
theory (plus signs), close to the melting temperature, for the particular value of the Lindcmann 
parameter c\ = 0.0153, for which the curve intersects the T-axis at Tm = 0.00701. The points lie 
on a straight line just as the prediction of the shaking theory, i.e., the dashed line Ft ^ 1/(T^ — T), 
yielding the value one for the critical exponent. 
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